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Bilateral Filter for Image Denoising 
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Abstract 

The bilateral filter has diverse applications in image processing, computer vision, and computational photog¬ 
raphy. In particular, this non-linear filter is quite effective in denoising images corrupted with additive Gaussian 
noise. The filter, however, is known to perform poorly at large noise levels. Several adaptations of the filter have 
been proposed in the literature to address this shortcoming, but often at an added computational cost. In this 
paper, we report a simple yet effective modification that improves the denoising performance of the bilateral filter 
at almost no additional cost. We provide visual and quantitative results on standard test images which show that 
this improvement is significant both visually and in terms of PSNR and SSIM (often as large as 5 dB). We also 
demonstrate how the proposed filtering can be implemented at reduced complexity by adapting a recent idea for 
fast bilateral filtering. 


Index Terms 

Image denoising, bilateral filter, box-filter, improvement, fast algorithm. 


I. Introduction 

Linear smoothing filters, such as the classical Gaussian filter, typically work well in applications where the 
amount of smoothing required is small. For example, they are very quite effective in removing small dosages of 
additive noise from images. However, when the noise floor is large and one is required to average more pixels 
to improve the signal-to-noise ratio, linear filters tend to to over-smooth sharp image features such as edges and 
corners. This over-smoothing can be alleviated using some form of data-driven (non-linear) diffusion, where the 
quantum of smoothing is controlled using the image features. A classical example in this regard is the anisotropic 
diffusion of Perona and Malik IT) . The insight of the authors was to take the standard diffusion equation and 
turn it into a non-linear differential equation by controlling its diffusivity using the gradient information. This 
automatically attenuated the blurring in the vicinity of edges. In practice, the associated differential equation is 
numerically solved using an iterative solver. While the Perona-Malik diffusion is known to be mathematically 
ill-posed 121, it is known to be numerical stable in practice and performs reasonably well on real data. A delicate 
aspect of this scheme is the choice of the stopping criteria which often critically determines the final result. 
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(a) Barbara. (b) Corrupted (cr = 20, 22.11 dB).(c) Standard Bilateral (26.62 dB). (d) Oracle Bilateral (30.05 dB). 

Fig. 1. Standard versus oracle bilateral filter. The standard and the oracle filters respectively use the comapted and the clean image to 
compute the range filter. Notice that the PSNR improves by more than 3 dB, and the image from the oracle looks visibly much better than 
that obtained using the standard filter. The result in (d) is the “best” one can hope to achieve if, instead of the clean image, some proxy 
is used to compute the range filter. 


A. Bilateral Filter 

The bilateral filter was proposed by Tomasi and Maduchi 13 as a simple, non-iterative alternative to the 
Perona-Malik diffusion. The origins of the filter can be traced back to the work of Lee H and Yaroslavsky 0. 
The SUSAN framework of Smith and Brady 0 is also based on a similar idea. 

For a discrete image (/(i))jgz2. the bilateral filter is given by 

,, Eien5a,(j) 5a,,(/(i-j) -/(i)) /(i-i) 

In this formula, 5crs(i) is the spatial filter defined on some neighbourhood ft and garit) is the range filter. 
Typically, U is a square neighbourhood, U = [—lU, lU] x [—IF, lU], and both the spatial and range filters are 
Gaussian; 

and (f) = exp 

The support W of the spatial filter is usually set to be Sus- 

The spatial filter puts larger weights on pixels that are close to the pixel of interest compared to distant pixels. 
On the other hand, the range filter operates on the intensity differences between the pixel of interest and its 
neighbors (which makes the overall filter non-linear). The role of the range filter is to restrict the averaging to 
neighbouring pixels whose intensities are similar to that of the pixel of interest. In particular, f{i — j) — f{j) in 
O is close to zero when both i and its neighbors {i — j : j G 0} belong to a homogenous region. In this case, 
— j) — /(i)) ~ 1, and ([T]i effectively acts as a standard Gaussian filter. On the other hand, consider the 
situation in which the pixel of interest i is in the vicinity of an edge. If i — j and i are on the opposite sides 
of the edge, then gar{f{t — j) — /(i)) relatively small compared to what it is when i — j and i are on the 
same side of the edge. This effectively prohibits the mixing of pixels from different sides of the edge during 
the averaging, and hence avoids the blurring that is otherwise induced by linear filters. 
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B. Present Contribution 

The bilateral filter has found widespread applications in image processing, computer vision, and computational 
photography. We refer the interested reader to Q and the references therein for an exhaustive account of various 
applications. 

Our present interest is in the image denoising applications of the filter IH, 13, ifTOl . The bilateral filter has 
received renewed attention in the image processing community in the context of image denoising ifTTH . am . It 
is well-known that, while the filter is quite effective in removing modest amounts of additive noise, its denoising 
performance is severely impaired at large noise levels Q, 113 . To overcome this drawback, different iterative 
forms of the filter were proposed in |[8], for example. In a different direction, it was shown by Buades et al. 
03 that a patch-based extension of the filter can be used to bring the denoising performance of the filter at 
par with state-of-the-art methods. However, this and other advanced patch-based methods iflTll . ITSll . IThll are 
much more computation-intensive than the bilateral filter. 

In this paper, we demonstrate how the denoising performance of the bilateral filter can be improved at 
almost no additional cost by incorporating a simple pre-processing step into the framework. To the best of 
our knowledge, this improvement has not been reported in the literature on bilateral filtering-based denoising. 
Although the present improvement is not of the order of the improvement provided by K-SVD and BM3D, 
we will demonstrate that the improved bilateral filter is often competitive with the Non-Local Means (NLM) 
filter lfT3 . while being significantly cheaper. Moreover, we also describe a fast algorithm for the modified filter 
which should be of interest in real-time denoising. 

C. Organization 

The rest of the paper is organized as follows. In Section [B] we briefly describe the denoising problem and 
the standard metrics that are used to quantify the denoising performance. We then introduce the proposed 
improvement of the bilateral filter in the context of denoising. In this section, we also report a fast algorithm 
for the proposed filtering. Experimental results on synthetic and natural images are provided in Section |nll and 
we conclude the paper in Section m 


II. Improved Bilateral Filter 

We consider the problem of denoising grayscale images that are corrupted with additive white Gaussian noise. 
In this setup, we are given the corrupted (or noisy) image 

/(i) = /o(i) + cr • ni (i G /), (2) 

where 

• / is some finite rectangular domain of 7?, 

• (/o(i))ie/ is the unknown clean image, and 

• (njig/ are independent and distributed as Af{0, 1). 

The goal is find a denoised estimate /(i) of the clean image from the corrupted samples. The denoised image 
should visually resemble the clean image. To quantify the resemblance, two standard metrics are widely used 
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in the image processing literature, namely the peak signal-to-noise ratio (PSNR) and the structural similarity 
index (SSIM) HT). The PSNR is defined to be 10logio(255VMSE), where 

' ' le/ 


A. Proposed Improvement 

In linear diffusion, the clean image is estimated by linearly averaging the noisy samples. The averaging 
process successfully brings down the noise floor in homogenous regions by a factor of s/L, where L is the 
length of the filter ifTSl . However, the filter also implicitly acts on the underlying clean image in the process. 
As a result, it introduces blurring in the image features besides reducing the noise floor. This can precisely be 
overcome by applying the bilateral filter on the corrupted image m, 13 , QO). In this regard, note that the range 
filter operates on the noisy samples. In other words, the corrupted image is used not just for the averaging 
but also to control the blurring via the range filter. What if the range filter could directly operate on the clean 
image? That is, instead of O, suppose we consider the formula 

Ejgng^.(j) g^r(/o(i-j) -/o(i)) /(i-j) 


/obf(i) = 


(3) 


Ejen 9<ts U) g<T,.(/o(i - j) - /oW) 

The denoising result obtained using this “oracle” filter is compared with that obtained using ([T]| in Figure [T] It 
is not surprising that the result obtained using the oracle filter is visibly better and has higher PSNR. 

Of course, the problem is that we do not have access to the clean image in practice, and thus the oracle 
bilateral filter cannot be realized. Nevertheless, one could consider some form of proxy for the clean image. 
For example, one could use an “iterated” bilateral filter |0 where the output of the bilateral filter is used as 
a proxy. However, this requires us to compute O twice, which doubles the run time of the filter. Our present 
proposal is simply to use the box-filtered version of the noisy image as a proxy. In other words, the proposed 
improved bilateral filter (in short, IBF) is given by 

Ejgn g-T. (g) g<T.(/(i - j) - /(i)) /(i - j) 


/ibf(i) = 


T,jen9^si3) g<T,.(/(i-g)-/W) 


where 


/(i) = 


(2L + 1)2 

^ ^ ^ je[-L,L] 


/(i-g)- 


(4) 


(5) 


Clearly, the amount of smoothing induced by the box-filter is controlled by L. When L is very small, /(i) Ri f{i), 
and (S behaves as O- At the other other end, the image structures are over-smoothed when L is large and 
this makes /(i) a bad proxy for the original image. Thus, L should not be too small and neither too large. We 
will report the appropriate choice of L in the sequel. 


B. Fast Implementation 

The cost of computing (|4]i is almost identical to that of computing ([T]l, since the additional cost of computing 
(|5]l is negligible in comparison. More precisely, the cost of computing (|4|i is 0(|cr^|) per pixel, since the support 
H of the spatial filter is proportional to cr^. On the other hand, it is well-known that the box-filter in Q can 
be computed using 0(1) operations per pixel lfT9l . We now explain how we can implement Q using 0(1) 
operations (with respect to as) as a straightforward extension of the algorithm proposed in ifTSl . Il20l . 
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Fig. 2. Approximation of the range kernel on the interval [—255,255] using raised-cosines of different orders. In this case, 

= 30 . Note that the approximation around the origin is already good when N = b, but the approximation at the boundaries is far 
from satisfactory (due to the periodic nature of the raised-cosine). Higher-order raised-cosines are required to suppress the oscillation at 
the boundaries and to force non-negativity and monotonicity. 


For completeness, we present the main ideas behind the fast algorithm in ifTSl . Note that the effective domain 
of the range filter in (|4|i is the interval [—T,T], where T is the maximum value of /(i — j) — /(i) over all 
1 G I and j G [—W, W]'^. In other words, T is the maximum “local” dynamic range of /(i) over square boxes 
of length 2W. Note that the complexity of computing T is comparable to that computing the filter, namely 
O(IV^) operations per pixel. A fast algorithm for computing T was however later proposed in ll20l . which we 
will use in this paper. 

It was observed in ifTSl that 


gaAt) = Jim 

N—¥oc 


cos 


t 


N 


(6) 


The functions on the right are called raised cosines, and we refer to N as its order. Note that while ® holds for 
every t, we only require a good approximation for t G [—T,T]. Moreover, the raised cosine should ideally be 
positive and monotonic on this interval. In particular, one can verify that if N is larger than Nq = 0.405(r jcTr)^, 
then the raised-cosine is positive and monotonic on [—T, T], In other words, any raised-cosine of order N > Nq 
is an acceptable approximation of the Gaussian range filter. The approximation process in demonstrated in Figure 

E 


Now, using the identity 2 cos 9 = 
expansion, we can write 




(where i denotes the imaginary unit ^ —V) and the binomial 


cos 


r^fNJ. 


n N 


N 


n—0 


where 


1 


= '^Cn exp {lUJnt) , 
(2n - N) 


Cn — 2 ^^ ^ ^ and u}ri — 


(7) 


( 8 ) 


ars/N 

Plugging (|7l) into (IHi, using the multiplication-addition property exp(u -f u) = exp(u) exp{u), and exchanging 
the summations, we can express the numerator of (lUi as 


N 


c„exp (-iw„/(l)) {Fn*gaJ{l), 

n—0 


(9) 
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where 

= exp(ia;„/(i))/(i), 

and 

{Fn * 5aJ(l) = X! 9aAj)Fn{y - j) 
jen 

denotes the Gaussian filtering of Fn{i). Similarly, we have the following approximation for the denominator of 

(Eli: 

N 

^ C„exp {-lUJnfil)) {Gn * (10) 

n—0 

where 

G„(i) = exp(za;„/(i)). 

To summarize then, by using the raised-cosine approximation of the Gaussian range filter, we can express 
the numerator and denominator of (El as a linear combinations of Gaussian filters applied on the images Fn (i) 
and Gra(i). It is well-known that Gaussian filtering can be computed using 0(1) operations per pixel (i.e., 
independent of as) using recursions ifTOll . As a result, we can compute (13l and (fTOl . and hence the overall filter 
(El, using 0(1) operations per pixel. 

In fact, it is further possible to cut down the run time without appreciably sacrificing the approximation 
using truncation ll20ll . In particular, it can be verified that the contribution of the central terms in (l2| to the 
overall approximation is less compared to the other terms. Indeed, the distribution of the normalized binomial 
coefficients (c„) is bell-shaped with a peak around N/2. As a result, one can truncate the sum away from the 
central peak and tradeoff speed versus accuracy. In particular, given some tolerance £ > 0, we incrementally 
find the largest M such that 

Cm + ... + cm-m >1 — e/2. (11) 

We can then further approximate (|7| using the truncated sum 

N-M 

c„ exp (iWnt). (12) 

n—M 

Note that the error between (fT^ and ([Tl) is 

N N-M 

Cn exp {lUJnt) - X! Cn exp (tUJnt) , 
n—0 n—M 

whose magnitude is, by construction, at most as large as 2(co -I-... -I- cm) < £■ Using (fT^ . we can further cut 
down the run time of (El by a factor of about 2M/N, without any appreciable change in denoising performance. 
For large N, one can test for the condition cm + ... + cn-m > 1 — s/2 without having to compute all the c„. 
In particular, the Chernoff-inequality for the binomial distribution gives us the estimate 

M = \ (W- v'4iVlog(2/£)) , (13) 

which is quite accurate when N > 100 ll20l . 

The complete algorithm is summarized in Algorithm [1] Here, we have used [N] to denote the set of numbers 
0,1,..., A^, and [N/M] to denote the set of numbers 0,..., M and N — M,... ^ N. In step 6 (c), we have 
used G* (i) to denote the complex conjugate of G„(i). 
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Data; Corrupted image /(i), and filter parameters as,crr,e. 

Result: Filtered image /ibf(i)- 

Initialize; Images P(i) = 0 and (5(i) = 0 for all i; 

Compute box-filtered image /(i) using 0; 

Compute local dynamic range T of /(i); 

N = 0.405(r/crr)2; 
switch do 

case < 40 
M = 0; 

Set Cn,uJn for n £ [N] using ®; 

endsw 

case 40 < < 100 

Set Cn,uJn for n £ [A^] using ®; 

Fix £ = 0.01 and find largest M such that (fTTll holds; 
endsw 

case A^ > 100 

Set M using (fOl) with £ = 0.1; 

Set Cn,uJn for n £ [N/M] using ®; 

endsw 

endsw 

for n = M,... ,N — M do 

0 ( 1 ) = exp(*a;„/(i)); 

P(i) = G(i)/(i); 

P(i) = c„G*(i); 

Q(i) = Q(i) + iF(i)-(G*p<,J(i); 

end 

/ibf(i) = P(i)/Q(i) for all 1 ; 

Algorithm 1: Fast Improved Bilateral Filtering 


III. Experiments 

A. Complexity and Run Time 

The complexity of the direct implementation of the proposed filter is identical to that of the standard bilateral 
filter, namely O(tTs). On the other hand, for an image with maximum local dynamic range T, the complexity 
of the fast implementation proposed in Section Hl-BI is 0 (T^/ct^). The final run time is however determined 
by the constants that are implicit in the above complexity estimates. In practice, the main speed-up is due to 
the fact that Gaussian filtering can be implemented very efficiently using standard packages. For example, the 
image filtering in step 6 (d) and (e) can be done using the optimized “imfilter” routine in Matlab. In Table |I] we 
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compare the run times of the fast implementation and the direct implementation (for typical filter parameters) 
on the Barbara image. 


TABLE I 

Comparison of the average run times of the direct and the fast implementation of the proposed eilter eor 

DIEEERENT (os, CTr). WE USED Barbara EOR THE COMPARISON AND THE NOISE LEVEL WAS SET AT CT = 20. FOR BOTH 
IMPLEMENTATIONS, WE TOOK THE SUPPORT OP THE SPATIAL GAUSSIAN TO BE VC = Scts . 


(2,15) 

(4,20) 

(3,25) 

(5, 30) 

(3, 35) 

(4, 40) 

Direct 16.5s 

60.5s 

35.3s 

93.8s 

35.5s 

60.5s 

Fast 0.52s 

0.61s 

0.52s 

0.47s 

0.43s 

0.47s 


All computations were performed using Matlab on a quad-core 2.70 GHz Intel machine with 16 GB memory. 
It is clear from the table that a significant acceleration is achieved using the fast algorithm. In particular, notice 
that the fast implementation takes about 0.5 seconds for different parameter settings. On the other hand, notice 
the run time of the direct implementation scales up quickly with the increase in the width of the spatial filter. 
We note that the run time of the fast implementation can be further cut down using a parallel (multithreaded) 
implementation of step 6 in Algorithm [1] 

B. Optimal Choice of L 

We now come to question about the choice of the optimal length L for the box-filter in (|5l). We performed 
exhaustive some simulations in this direction, the results of some of which are reported in Figure [3 For these 
simulations, we conclude that a box-filter with L = 1 (3 x 3 blur) is optimal for most settings. A possible 
way to explain this is that this small filter is able to suppress the noise without excessively blurring the image 
features. 

C. Denoising Results 

We now present some denoising results to demonstrate the superior denoising performance of the proposed 
filter. In Figure 0 ] we compare the proposed filter with the standard and the oracle filter on a synthetic test 
image. Notice that the improvement in PSNR over the standard filter is more than 10 dB. This does not come 
as a surprise since this particular test image has a lot of sharp intensity transitions. While the bilateral filter is 
already known to work well for such images, what this result demonstrates is that we can further improve its 
performance using the proposed modification. Moreover, note the PSNR of the proposed filter is close to that 
of the oracle bilateral filter (which uses the clean image to compute the range filter). In Figure |5] we compare 
the proposed filter with the iterated bilateral filter. Notice that while the PSNR from the iterated filter is within 
a dB of that obtained using the proposed filter, the denoised image from the latter looks visibly better. This is 
also confirmed by the SSIM indices reported in the figure. 

We next compare the denoising performance of the proposed modification with the standard bilateral filter on 
certain standard test images ||2TI| . The PSNR and SSIM indices of the proposed modification and the standard 
bilateral filter are reported in Table In] We independently optimize the standard and the improved bilateral filter 
with respect to CTs and (Jr- Also, we use L = 1 for the improved bilateral filter. Notice that the proposed filter 
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(a) Barbara. (b) Lena. (c) House. 


Fig. 3. The plots show how the PSNR varies with L at various noise levels. For all the experiments, we used cjs = 6 and <Jr = 20. 
Notice that the PSNR is consistently highest when L = 1. 


starts to perform better beyond a certain noise level (cr Ri 20). On the other hand, the SSIM improvement is 
already evident for all the images beyond cr = 15. This is because, at low noise levels, the proposed box-filtering 
does more blurring than noise suppression, which brings down the overall signal-to-noise ratio. Indeed, when 
the noise floor is small, the corrupted image is already a good proxy for the clean image. However, notice that 
the improvement in SSIM is quite significant for all the images at large noise levels, and the PSNR improvement 
is often as large as 5 dB. For a visual comparison, some of the results of the denoising experiments are shown 
in Figures |6l |2l and [H In Table [III we compare the proposed filter with some sophisticated image denoising 
methods cited in the introduction, namely, NLM lfT3l . K-SVD ITSl . and BM3D IfTO . NLM is essentially a 
patch-based extension of the bilateral filter, in which image patches (groups of neighbouring pixels) are used 
for comparing neighbouring pixels. The latter methods are based on sparse-coding and collaborative filtering 
and are significantly more sophisticated. What we find quite interesting is that, beyond a certain noise level, 
the proposed filter is competitive with NLM for most of the test images in Table [Hi except for Barbara and 
Cameraman. The denoising performance is in general inferior to K-SVD and BM3D. We do not find this 
surprising since (among other things) they heavily rely on the sparsity of natural images (in appropriate bases) 
to improve the PSNR by few extra dBs. However, we believe that the present work is relevant in the context of 
the recent results in and lEl, where the bilateral filter is used to obtain state-of-the-art denoising results. 

IV. Conclusion 

We demonstrated that by using a box-filtered image in the range filter, one can substantially improve the 
denoising performance of the bilateral filter, at almost no additional cost. While the basic idea is quite simple, 
it is nevertheless quite effective in improving the denoising performance of the filter. Exhaustive denoising 
results on test images were provided in this direction. This address a well-known pathology of the bilateral 
filter, namely, that its denoising performance begins to degrade quickly with the increase in noise level. An 
interesting finding was that the proposed filter is often competitive with the computation-intensive non-local 
means filter. We also presented a fast algorithm for the proposed filter that can dramatically cut down the run 
time. As future work, we plan to investigate how one can combine the standard and the proposed filter so as 
to consistently obtain the best denoising performance at all noise levels. 
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TABLE II 

Comparison of the Standard Bilateral Filter (SBF) and the Improved Bilateral Filter (IBF) in terms of PSNR and 
SSIM AT NOISE LEVELS a = 10, 15, 20, 25, ... , 55, 60. WE ALSO COMPARE THE PSNR WITH THOSE OBTAINED USING MORE 
SOPHISTICATED METHODS SUCH AS NON-LOCAL MEANS (NLM) 1731 . K-SVD 1731 . AND BM3D (76). FOR A FIXED IMAGE AND 
NOISE LEVEL, WE TUNE Us AND Ur TO INDEPENDENTLY OPTIMIZE THE PSNRS OBTAINED USING SBF AND IBF. WE ALSO TUNE THE 
PARAMETERS OF NLM TO GET THE OPTIMAL PSNR AT EACH NOISE LEVEL. IF THE PSNR OBTAINED USING IBF IS HIGHER THAN 
THAT OBTAINED USING NLM, WE MARK THE EORMER IN BOLDFACE. 


Image 

Filter 





PSNR (dB) 






SBF 

31.41 

28.78 

27.13 

25.51 

23.63 

21.93 

20.30 

18.76 

17.40 

16.21 

15.10 


IBF 

25.60 

25.45 

25.18 

24.85 

24.47 

24.15 

23.89 

23.60 

23.42 

23.24 

22.96 

Barbara 

NLM 

33.04 

30.86 

29.27 

27.99 

27.16 

26.42 

25.74 

25.04 

24.50 

24.04 

23.67 


K-SVD 

34.42 

32.27 

30.76 

29.44 

28.40 

27.43 

26.61 

25.87 

25.23 

24.72 

24.18 


BM3D 

34.98 

33.05 

31.68 

30.55 

29.63 

28.76 

27.88 

27.63 

26.99 

26.54 

26.13 


SBF 

33.58 

31.60 

29.74 

27.30 

24.78 

22.70 

20.80 

19.15 

17.66 

16.36 

15.25 


IBF 

33.27 

32.49 

31.49 

30.59 

29.79 

29.19 

28.66 

28.05 

27.62 

27.17 

26.62 

Lena 

NLM 

34.06 

32.33 

30.99 

29.89 

29.07 

28.39 

27.74 

27.18 

26.72 

26.23 

25.85 


K-SVD 

35.50 

33.65 

32.40 

31.27 

30.44 

29.67 

29.01 

28.38 

27.79 

27.31 

26.90 


BM3D 

35.88 

34.20 

33.02 

32.03 

31.16 

30.47 

29.76 

29.44 

28.96 

28.57 

28.16 


SBF 

33.73 

33.70 

29.64 

27.18 

24.67 

22.62 

20.74 

19.06 

17.61 

16.36 

15.22 


IBF 

33.17 

32.39 

31.38 

30.50 

29.79 

29.11 

28.49 

27.72 

27.26 

26.82 

26.28 

House 

NLM 

34.63 

33.00 

31.63 

30.56 

29.34 

28.47 

27.69 

26.93 

26.36 

25.70 

25.00 


K-SVD 

35.89 

34.39 

33.07 

32.15 

31.29 

30.36 

29.59 

28.89 

27.90 

27.28 

27.08 


BM3D 

36.80 

35.05 

33.71 

32.86 

32.15 

31.30 

30.88 

30.16 

29.99 

29.12 

28.69 


SBF 

32.97 

30.74 

28.80 

26.54 

24.14 

22.01 

20.20 

18.64 

17.32 

16.08 

15.00 


IBF 

31.31 

30.61 

29.80 

28.72 

27.93 

27.03 

26.34 

25.69 

25.14 

24.65 

24.22 

Peppers 

NLM 

32.91 

30.71 

29.23 

28.06 

27.14 

26.26 

25.51 

24.98 

24.40 

23.82 

23.35 


K-SVD 

34.27 

32.34 

30.87 

29.73 

28.87 

28.10 

27.33 

26.79 

26.13 

25.66 

25.00 


BM3D 

34.70 

32.68 

31.27 

30.21 

29.21 

28.54 

27.67 

27.23 

26.75 

26.26 

25.89 


SBF 

32.03 

29.90 

28.40 

26.39 

24.21 

22.24 

20.52 

18.93 

17.50 

16.26 

15.17 


IBF 

29.96 

29.55 

28.90 

28.15 

27.46 

26.90 

26.45 

25.93 

25.56 

25.16 

24.81 

Boat 

NLM 

31.93 

29.93 

28.57 

27.60 

26.90 

26.25 

25.68 

25.12 

24.58 

24.19 

23.84 


K-SVD 

33.63 

31.69 

30.36 

29.24 

28.42 

27.67 

27.01 

26.46 

25.95 

25.49 

25.07 


BM3D 

33.88 

32.09 

30.79 

29.81 

29.06 

28.35 

27.69 

27.10 

26.73 

26.27 

25.98 


SBF 

32.65 

30.25 

28.54 

26.33 

24.23 

22.24 

20.44 

18.91 

17.40 

16.19 

15.08 


IBF 

27.58 

27.35 

26.98 

26.49 

25.91 

25.41 

24.96 

24.63 

24.27 

23.89 

23.62 

Cameraman 

NLM 

32.61 

30.00 

28.57 

27.70 

27.01 

26.26 

25.39 

24.75 

24.28 

23.92 

23.22 


K-SVD 

33.73 

31.37 

29.96 

28.87 

28.00 

27.33 

26.69 

26.30 

25.74 

25.16 

24.92 


BM3D 

34.16 

31.84 

30.41 

29.53 

28.60 

27.84 

27.07 

26.63 

25.98 

25.72 

25.38 

SSIM (%) 

Barbara 

SBF 

89.51 

82.24 

75.21 

64.01 

51.54 

42.2 

35.15 

29.51 

24.77 

21.29 

18.39 

IBF 

76.36 

75.16 

73.00 

69.87 

67.61 

65.44 

63.71 

61.97 

59.87 

59.30 

56.51 

Lena 

SBF 

88.45 

84.12 

74.56 

60.29 

45.20 

34.22 

26.54 

21.18 

16.69 

13.80 

11.55 

IBF 

87.80 

86.27 

83.52 

81.68 

79.89 

78.51 

77.14 

75.38 

72.84 

71.49 

68.01 

House 

SBF 

86.92 

83.50 

74.02 

59.27 

45.17 

34.45 

27.19 

21.78 

17.69 

15.04 

12.86 

IBF 

86.45 

84.94 

82.37 

81.33 

80.12 

78.78 

76.89 

73.50 

72.19 

70.69 

67.83 

Peppers 

SBF 

90.86 

86.06 

78.37 

66.02 

52.56 

42.49 

34.40 

27.26 

23.14 

19.77 

17.23 

IBF 

89.83 

88.44 

86.01 

82.09 

80.29 

75.76 

73.52 

71.98 

72.03 

70.29 

66.98 

Boat 

SBF 

85.62 

79.31 

73.39 

62.01 

50.20 

38.88 

31.67 

25.85 

21.42 

17.97 

15.42 

IBF 

80.51 

79.47 

77.30 

74.20 

71.86 

69.30 

67.72 

65.61 

63.22 

61.70 

59.99 

Cameraman 

SBF 

88.91 

82.21 

75.15 

61.39 

49.36 

39.14 

32.06 

27.23 

23.05 

20.21 

17.59 

IBF 

84.07 

82.83 

80.09 

76.18 

73.76 

73.83 

73.61 

68.85 

69.46 

64.46 

62.49 
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(a) Comapted (cr = 40,16.12 dB).(b) Standard Bilateral (18.92 dB). (c) Proposed (29.23 dB). (d) Oracle Bilateral (32.16 dB). 


Fig. 4. Denoising results for Checker (original size 150 x 150). Parameters used: cts = 3, ct^ = 30, £ = 0.01, and W = 3f7s. Notice 
that the image in (c) looks significantly better than that in (b) and is similar to the oracle in (d). 



(a) Con-upted (a = 25, 20.13 dB).(b) Standard Bilateral (25.52 dB). (c) Proposed (29.92 dB). (d) Iterated Bilateral (28.71 dB). 

Fig. 5. Denoising results for House. The filter parameters in this experiment are identical to the ones reported in Figure[4] For the iterated 
bilateral filter, we use (b) to compute the range filter instead of using (a), however the averaging is performed on (a). Thus the run time 
of the iterated filter is roughly twice that of the proposed filter. Notice that (d) has some granular artefacts particularly along the edges, 
even though its PSNR is within a dB of that of (c). In this regard, we note that the SSIM index of (c) is 0.82, while that of (d) is 0.76 



(a) Lena {a = 30,18.57 dB). (b) SBF (24.72 dB, = 2, cr^ = 40). (c) IBF (29.82 dB, Cs = 2,(Tr = 20). 

Fig. 6. Denoising results using standard bilateral filter (SBF) and improved bilateral filter (IBF). We tuned the parameters of SBF and 
IBF to get the optimal PSNR in either case. The parameters settings are indicated in the respective captions. 
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(a) Peppers (a = 30,18.55 dB). (b) SBF (24.18 dB, cts = 2,ar = 40). (c) IBF (27.91 dB, cr^ = 2,ar = 20). 

Fig. 7. Denoising results using standard bilateral filter (SBF) and improved bilateral filter (IBF). We tuned the parameters of SBF and 
IBF to get the optimal PSNR in either case. The parameters settings are indicated in the respective captions. 



(a) Cameraman (cr = 35,17.24 dB). (b) SBF (22.30 dB, cr^ = 2, cr^ = 40). (c) IBF (25.45 dB, cr^ = 2, cr^- = 20). 

Fig. 8. Denoising results using standard bilateral filter (SBF) and improved bilateral filter (IBF). We tuned the parameters of SBF and 
IBF to get the optimal PSNR in either case. The parameters settings are indicated in the respective captions. 
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